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Abstract 

In this article, we retain the basic idea and at the same time generalize 
Kawasaki's dynamics, spin-pair exchange mechanism, to spin-pair redistribu- 
tion mechanism, and present a normalized redistribution probability. This 
serves to unite various order-parameter-conserved processes in microscopic, 
place them under the control of a universal mechanism and provide the ba- 
sis for further treatment. As an example of the applications, we treated the 
kinetic Gaussian model and obtained exact diffusion equation. We observed 
critical slowing down near the critical point and found that, the critical dy- 
namic exponent z = 1/v = 2 is independent of space dimensionality and the 

assumed mechanism, whether Glauber-type or Kawasaki-type. 
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Irreversible dynamic systems exhibit complicated but interesting nonequilibrium phe- 
nomena near the critical point. In spite of their complexity, the interesting dynamic critical 
behaviors have been attracting a lot of researchers for many decades. Within the vast 
body of literature, two pioneering works completed by Glauber |1| and Kawasaki g] have 
been regarded as a milestone. Great progress has been achieved with the application of 
Glauber's single-spin flip mechanism and kawasaki's spin-pair exchange mechanism. These 
two schemes catch the inherent essential process bound in. 

Kawasaki's dynamics deals with a system consisting of an array of iV coupled spins. The 
coupling between spins is represented by a set of transition probabilities of spin exchange. 
Initially, the focus was on Ising model. The spins may exchange with their nearest neighbors, 
and in this way, the system evolves while the total spin remains conserved. In later studies 
the idea of exchange has been proved very successful as it catches the essential nature. It has 
important applications in Ising model |3],f|] and Ising-like models such as lattice gas model 
0-0, Blume- Emery-Griffiths (BEG) model ||[!| and others |I0| , |r[j| . These applications have 
contributed a lot to the understanding of the thermodynamics, and have erected itself as 
the basic mechanism in the class of order-parameter-conserved processes. 

However, its initial embodiments was closely tied up with the simpleness of Ising model, 
in which the spins can only take two values, ±1, and there exists only nearest-neighbor 
interactions. It turned to be limited when applied to systems other than Ising or Ising-like 
models. The same situation exists in Glauber's dynamics. Recently Zhu and Yang success- 
fully generalized Glauber's single-spin flip mechanism to single-spin transition mechanism 
and gave a normalized version of the transition probability |12[]. The applications in sev- 



eral models have yielded encouraging results |T2|JT3|| . In this article, along the same line, 
we retain the basic idea and at the same time generalize Kawasaki's dynamics, spin-pair 
exchange mechanism, to spin-pair redistribution mechanism, and also present a normalized 
redistribution probability. 

The limit is to be removed and the spins may take various discrete values (discrete-spin 
model) or continuous values (continuous-spin model). Because in Ising model spins can only 



have two opposite values, simply direct exchange may be enough to describe the way the 
system evolves. In other different models, such as Gaussian model, Potts model, or XY 
model, this simple picture may not be as capable. On the other hand, as is mentioned 
above, the conservation of the order-parameter has been regarded as the most important 
feature of this class of processes (and a necessary result of the exchange mechanism). Based 
on these considerations, in the generalized mechanism, spin pair redistribution mechanism, 
two neighboring spins, aj., cjh-Ij 110 longer merely exchange with each other. Instead, they 
may take any values, a^, cfk+i, as long as their sum remains conserved, or, their sum is 
redistributed; and we call this, spin-pair redistribution. 

The probability distribution function P (01, . . . , a^] t), or simply P ({a} ; t), denotes the 
probability of the TV-spin system being in the state (o"i, . . . , <tn), or simply {a}, at time t. 
The jth and Zth spins are two neighboring spins, and let Wji (<Jj<Ji — > be the probability 
per unit time that they are redistributed while the others remain unchanged. Then, on the 
supposition of neighboring spin-pair redistribution, we may write the time derivative of the 
function P ({a} ; t) as 

jP (M ;*)=EE {-W fl - &&) P ({a} ; t) 

(Jl) &jA 

+ Wji -> OjOi) P ({(Ti^y} , frj, au t)} (1) 

This is a probability equation, in which the first term in the right-hand side of (|l|) denotes 
the decrease of the probability distribution function P ({a} ; t) per unit time due to the 
redistribution of the spin pair from initially OjO\ to various values OjO^, and the second term 
denotes the increase of P {{&} ; t) due to the redistribution of the spin pair from initially 
various values OjO\ to finally cjjCT;. (Clearly, in Ising and Ising-like systems, it does become 
Kawasaki's picture.) We shall refer to ([I]) as the master equation since its solution would 
contain the most complete description of the system available. 

The key step for solving the master equation is the determination of the redistribution (or 
exchange, flip, transition, etc.) probability. Usually it can not be uniquely determined by the 
detailed balance condition, and thus, in Kawasaki's as well as Glauber's pioneering works, 
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some arbitrariness remained. We hope to make our choice of the spin-pair redistribution 
probability able to contain Kawasaki's in the specific Ising model, and applicable to various 
other spin systems, while at the same time appropriate, clearer, and more definite. Now we 
consider it in both mathematical and physical aspects. In mathematics, generally speaking, 
the probability must be positive and able to be normalized; in physics, we often require 
that a system in thermodynamic equilibrium satisfy the detailed balance condition; we 
also require the probability be ergodic as long as the total spin keeps conserved through 
the redistribution. Based on these considerations, for arbitrary neighboring jth and Ith 
spins, we can choose the redistribution probability Wji [pj^i — ► Vjvi) to satisfy the following 
conditions: 

(a) Ergodicity, positivity, and conservation of spin: 

Wji ((JjOi -> GjGi) > 0, Vfrj + a i = Gj + Gi 
Wji (GjGi -> GjGi) = 0, V(7j + 07 7^ Gj + Gi 

this may suggest a 5 function in the expression; 

(b) Normalization: 

(c) Detailed balance: 

Wji {GjGi — » GjGj) _ P eq (ffj, . ■ ■ , Gj, Gj, . . . , a N ) 
Wji (GjGi — > GjGi) P eq (<7i, . . . , Gj, Gi, . . . , G N ) ' 

in which 

P eg = \ exp \-m ({a})] , Z = ]T exp ({a})] , 

where P eg is the equilibrium Boltzmann distribution function, Z is the partition function, 
and Ti ({g}) is the system Hamiltonian. 

Although the spin pair redistribution probabilities still are not determined uniquely by 
the above restriction, there is less room left. The consideration that we have is similar to 
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that in the generalization of Glauber's dynamics ||12|| , that the redistribution of a neighboring 
pair depends merely on the momentary values of the spins surrounding them, or, the spins 
of the system, and the influence of the heat bath. Based on this, we can similarly assume 
that the redistribution probability Wji {cfjCfi — > tijd'i) only depends on the heat Boltzmann 
factor of the system. 

Wjl {cTjUl — > CTjCTl) OC 5 aj 



or 



Wji (ajai ajai) = —5 aj>avA . A exp 



-f3Hji (aj, a h {a m } 



(2) 



Using the normalization condition (b), we can write the normalization factor Qji as 



Qjl — ^aj,(Ti;aj,ai GXp 



In fact, the Hamiltonian coming from the interaction of spins unrelated to <jj, <Ji will be 
cancelled, and thus in the actual calculation, one needs only write the effective Hamiltonian. 
Obviously, Qji is related to the temperature, the surrounding spins and the conserved sum 

Of (Tj, 07. 

Compared with Kawasaki's expression, (fj) is a normalized version. In his expression, 
there is an a appearing in the exchange probability, which was assumed to be a constant. In 
Ising model, redistribution is in fact exchange and actually our expression is only a definite 
selection for constant a by extra restriction and physical considerations. 

Usually, we are interested in local magnetization. It is defined as 



q k (t) = (<r k (*)) = $> fc P ({a} ; t) . 



(3) 



According to the definition ([3D and the master equation ([[]), and using the normalization 
condition (b), the time-evolving equation of (t) can be derived as (see Appendix |A]) 

d 



dt 



q k (t) = -2dq k (*) + ££ 



P({a};t) 



(4) 



where d is the dimensionality of the system, and J2 W means summation over the nearest 
neighbors (clearly it is related to the dimensionality, too). 

Kawasaki's exchange mechanism was initially designed for the study of diffusion constant, 
and he himself obtained an approximate result for Ising model by first deriving an expression 
of spin flux [Q. Now we will treat the kinetic Gaussian model as an application of our 
redistribution mechanism, while our method is a direct one. 

The Gaussian model, proposed by T. H. Berlin and M. Kac, at first in order to make an 
Ising model more tractable, is an continuous-spin model. It has the same Hamiltonian form 
as the Ising model (three dimensional), 

N 

-(3H = K £ £ 

i,j,k=l w 

where J2w means summation over nearest neighbors. Comparing it with the Ising model, 
there are two extensions: First, the spins can take any real value between (—00, +00). 
Second, to prevent the spins from tending to infinity, the probability of finding a given spin 
between er^ and + da^ is assumed to be the Gaussian-type distribution 

/ (a ljk ) da ijk = \/^ ex P (^-^ijk^j dvijk, (6) 

where b is a distribution constant independent of temperature. Although it is an extension 
of the Ising model, the Gaussian model is quite different. In the equilibrium case, on 
translational invariant lattices the Gaussian model was exactly solvable, and later as a 
starting point to study the unsolvable models it was also investigated with mean field theory 
and the momentum-space renormalizat ion-group method. 

In the 3D kinetic Gaussian model, the system Hamiltonian and the spin distribution 
probability are (|5|) and (|6|), respectively. In this case there are six combined terms in the 
3D type time-evolving equation of the local magnetization, Eqs.(|j), and the details of these 
complex calculations are in Appendix [FJ. Here we give only the results, 

E ^ijk^i,j,k;i+l,j,k {&ijk&i+l,j,k ~^ &ijk&i+l,j,k) 
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2(b + K) ^ ai ~ l j' k ai d-i,k + a i,j,k-l + a i,j,k+l + a ijk + a i+l,j,k 

~ (J i+l,j+l,k — &i+2,j,k — &i+lj-l,k ~ ^i+lj.fe+l — Cj+l.^fc-l) + b (<7ijk + Ci+lj,fc)] 



^ &ijkWij j k;i-l,j,k{ a 'ijk a 'i-l,j,k ~~ > &ijk&i~l,j,k) 

®ijk I*-*!— i-,j,k 

Z (0 + K) 



&ijk t& i,j-\-l,k 

^—^ [K (a iijM1 + a^ k + + + <x m ,,, fc + a ijk + a,,.,, 

~ a i,j+l,k+l — 0~i,j+2,k — °"i,j+l,fc-l — &i+lj+l,k ~ a i-l,j+l,k) + b (<7ijk + ^ij+l,*)] 



1 



Vijkt&i,]— l,h 



, , , , (fijjfc+l + &i,j+l,k + &i,j,k-l + &i-l,j,k + &i+l,j,k + Cij'fe + 0"t,j-l,fc 

Z (0 + A J 

— C"i,j-l,fe+l — &i,j-2,k ~ a i,j-l,k-\ ~ a i+lj-X,k ~ a i-l,j-X,k) + & (°ijfc + °"i,j-l,fc)] 



3vi 3 k+Vij,k+i,Vijk+Vij,k+i&ijkWij t k]iJ,k+l (&ijk<J'i,j,k+l —> 6'ijk6'i,j,k+l / 

1 



[K (<Ji+lj,k + &i,j,k-l + &i-l,j,k + <Ti,j-l,k + CiJ+l.fc + Cjjfc + ^ij.fc+l 



2 (6 + AT) 

~~ ^i+lj'^+l — &i,j,k+2 — — "ij+l.fc+l ~~ OiJ-l.fc+l) + & (°"ijfc + 0"ij,fc+l)j j 



^ ^^ijk+^i,j,k-i,^ijk+^i,j,k-i^ijk^i,j,k;i,j,k~l k &i,j,k-l y &ijk&i,j,k-l) 

^ijk )Ci,j,k—l 

= [K (*i + l, j>k + + tTi-lJ* + tnj-ij, + *i, j+ l,k + *i jk + ^k-l 

~< 7 i+l,j,k-l ~ cr i,j,k-2 — — &i,j+l,k-l ~ °"i,i-l,fc-l) + b (<?ijk + <Ti,j,k-l)] ■ 

Substituting them into the time-evolving equation of the local magnetization (H), we get 
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j t q ijk (t) 



1 



~~ ~~ (Sijfc — Qi,j,k-l)\ 



2{b + K) 

+ [(Qi+i,j,k - qijk) - (qijk - qi-i,j,k)\ + [(Qi,j+i,k - qijk) - (qijk - ?i,j-i,fc)]} 
K 

+ 2 (b + K) ^ ^ qi ~ 1 ' j ' k ~ ~~ 9i-i,j-i,fc) + (2<7i-i,j,fc - - 9i-2,j,fc) 

+2 (2q i+1 j jk - q i+ ij + i t k - <?j+i,j-i,fc) + (2c/i+ij,fc - ^jfc - 9i+2,j,fc) 

+2 (2qij_i k — q i: j_ lik+ i — + (2q i j_ ljk — qij k — q i j_ 2 ,k) 

+2 (2qij+i, k - qij+i^+i ~ qi,j+i,k-i) + (2qi,j+i, k - (fofc - q%j+2,k) 

+2 (2qij tk -i — ^i— i — gj+i^fc-i) + (2q i: j tk -i — qij k — q^j^-2) 

+2 (2qij :k+ i - qi+i,j, k +i - qi~ij,k+i) + ^c/^fc+i - - 9i,j,fc+2)] • (7) 

With lattice constant a we can transform the above equation to be 

j?(o = 2(6^) 6 ( v * +v ; +v ')9(*) 

;(2Vj + V 2 ) + (2VJ + V 2 ) + (2V 2 + Vj) 



It is of the form of a diffusion equation 

dq(t) 
dt 



2(b + K)' 

+ (2V 2 Z + Vj) + (2V 2 + V 2 ) + (2V 2 + V 2 )] g (f) 
--K )V 2 q(t). 



DV 2 q (t) , 



where 



(8) 



(9) 



b + K\6 

With the same treatment, we can easily obtain the diffusion equation of one and two 
dimensional kinetic Gaussian model. Here we give only the results, 



d 



da 2 



dt q ° {t) b+ K 12,/ 



771- K) V 2 q, 



(10) 



where d is the system dimensionality. As we know that the critical point for <i-dimensional 
Gaussian model is K c = J/k,BT c = b/2d, thus, the diffusion equation we obtained reveals 



that the diffusion will get much slower when it is near the critical point. The linear equations 
we obtained for a single spin can be directly solved, but the solution of the diffusion equations 
may already give us satisfying information. For example, in ID case, it is 

1 



q(x,t) 



g(£,0)e d£. 



2\ADvrt J-oo 

In a specific example of the diffusion of a Gaussian type packet, <? (£, 0) = e~^ 2 , one will 
obtain 



q(x,t) 



1 



1+t/r 



cxp 



X 



where 



T 



1+t/r 



n -l 



(11) 



1 


' Aa 2 






(HI 


AD ~ 


b + K 





(12) 



is the relaxation time. When K — > K c = b/2, D — > and r — > oo, and this is a typical 
critical slowing down phenomenon. With the correlation length critical exponent v — 1/2 
and the following dynamical scaling hypotheses 

e~ ir-T c rv~e, 

one can obtain the dynamic critical exponent z = 2. The same result can be obtained for 



2D and 3D model. In earlier study [T2|] the same result z — 1/v — 2 for any dimension- 
ality was obtained with Glauber-type mechanism. Thus we find that, in Kinetic Gaussian 
model, the critical dynamic exponent is independent of space dimensionality and the dynamic 
mechanism. 

Before conclusion, we may note that actually, although a systematic formulation has been 
lacked, efforts towards wider application of Kawasaki's dynamics have never been stopped. 
For example, in a Monte Carlo simulation of the three-dimensional ferromagnetic Heisenberg 
model |TJ]], Zhang has suggested that two neighboring spins, a, and o~j, may rotate with their 
conserved sum being the axis. He found that the new scheme enabled the system to evolve 
to thermodynamic equilibrium faster, and commented that it might be more favorable in 
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reality. This is just an successful exploration of the spin-pair redistribution mechanism, and 
there are many other such examples in non-equilibrium statistics, though assuming different 
forms. 

To summarize, in this article, we presented a systematic formulation of the Kawasaki- 
type dynamics: spin-pair redistribution. As a natural generalization of Kawasaki's exchange 
mechanism, the new dynamics gives the system more freedom while keeping the order- 
parameter conserved. We have analyzed the redistribution in physics, given the master 
equation and a normalized redistribution probability determined by the heat Boltzmann 
factor. The presentation of this probability, being the key of the whole formulation, makes 
the mechanism mathematically well-organized and physically meaningful. As mentioned 
above, in many works focused on specific lattice systems there are already ideas of "redis- 
tribution" , and these efforts often turned to be rather fruitful. However without a universal 
theoretical foundation, these works still remained tentative to some degree. This article com- 
pactly and systematically provides this foundation, upon which the generalized Kawasaki 
dynamics becomes universal and can be directly applied to microscopic systems. Without 
any extra requirements, it has its advantage compared with some earlier approaches covering 
the same ground, such as numerical Ginzberg Landau approaches. 

The formulation is compact in mathematics, while on the other hand it is also quite open. 
People are able to introduce other elements into it. For example, for decades there have 
been attempts of the application of a combined mechanism, Glauber type and Kawasaki type 
both with a probability ||15|| . We can easily give the mathematical form of the combined 
mechanism, single-spin transition plus spin-pair redistribution. In fact for Gaussian model 
we have already obtained exact and physically clear results, which will be reported in later 
paper as further application. For another example we may note the recent interest in the 
small-world network JIB] effect on transportation. One can also modify the redistribution 
mechanism and directly study these effects in non-equilibrium statistical dynamics. 

The formulation of the mechanism is the chief purpose of this article. As an example 
of its applications we obtained exact diffusion equations in kinetic Gaussian model, and a 
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temperature-dependent diffusion coefficient which becomes zero at the critical point. We 
observed critical slowing down phenomenon in diffusion process, and found that at least in 
this specific case the critical dynamic exponent is independent of space dimensionality and 
the dynamic mechanism. 

Up to this point we have successfully generalized Glauber's single-spin flipping mech- 
anism to single-spin transition mechanism, and Kawasaki's spin-pair exchange mechanism 
to spin-pair redistribution mechanism. These two generalizations are of similar mathemat- 
ical form and become counterparts of each other in nonconserved dynamics and conserved 
dynamics respectively. 

The diffusion in Gaussian model studied in this article is only an example of the problems 
previously we were unable to treat. Based on the redistribution mechanism (or transition), 
one can also easily write the master equation and redistribution probability for Potts model, 
XY model, Heisenberg model and many other types (in principle arbitrary), and obtain the 
evolution of local magnetization (and correlation function, etc.) in these systems. Though 
exact treatment may be difficult, other methods can be later applied based on this founda- 
tion. At the same time this mechanism can be almost directly applied to conserved processes 
other than those in spin-lattice models. For example, it can be used to study the relaxation 
of granular material under shaking (mass conserved), or the activities of particles in space 
(particle number conserved). A feasible way may be like this: first we divide the system 
into small units and define the parameter (and its range) and the system Hamiltonian. The 
redistribution is not necessarily limited to neighboring units (we can modify the master 
equation accordingly) and the probability may assume a form other than the heat Boltz- 
mann factor (but it has to be normalized). Then we write the master equation and the 
evolving equations of the parameters we are interested in. This serves to unite the various 
conserved processes in microscopic, place them under the control of a universal mechanism, 
and provide the basis for further treatment, either exact, approximation, or Monte Carlo. 
The same is for the single-spin transition mechanism. We think it is very important for the 
general progress in non-equilibrium statistical dynamics. 
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APPENDIX A: EQUATION DERIVATION 



d 
dt 



E E E [-^kWji (<Tj<Ti -> P ({a} ; t) 
+or k Wji (aj&i -> <Tj<Ti) P ({cr^,J , <7j, 07; i)] 



Taking summation over all {a}, one will find the terms of those pairs unrelated to a k cancel 
with each other. So one will only have 

d 



dt 



EE E [- a k(W kyk+w (a k a k+w ^> a k a k+w )) P {{a} -,t) 

{a} w (J k ,a k+W 

+c k W k:k+w (a k a k+w — > cr k cr k+w ) P (cti, . . . ,a k , a k+w , . . . , a N ; t)\ 



= -2dq k W + EE 
M «° 

= -2% (t) 



E a kW k ,k+w {^k^k+w — ► O'kO'k+w) P ■ ■ ■ ,&k, Vk+w, ■ ■ ■ , &N', t) 



+ E E &kW k>k+w (a k a k+w — > a k a k+w ) P (<7i, . . . , (X fc , <7fc +U) , . . . , <Xat; t) 

W (Tl,...,(T fc ,(T fc + m ,(T fc ,(T fc + m ,...,(T J v 



=2^)+ee 



E ^W^fc^+u) (Ck^k+w — > O'kO'k+w) 



P({a};t). 



APPENDIX B: CALCULATIONAL DETAILS 



We give the details for one of them, and the remaining terms follow the same way. The 
spin pair redistribution probability can be expressed as 
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^Vi,j,k;i+l,j,k (&ijk&i+lj,k ~^ &ijk&i+l,j,k) 
1 



-da 



a ijk+ a i + l,j,k^ijk+ a i + l,j,k 



Qi,j,k;i+l,j,k 

X exp {K l&ijk + &i t j+i,k + &i,j-l,k + a i,j,k-\ + &i,j,k+l) + &ijk&i+l,j,k 

+&i+l,j,k (&i+2,j,k + &i+l,j+l,k + 0"i+l,i-l,fc + &i+l,j,k+l + CTi+l,i,fc-l )]}, 

and there are similar expressions for Wij^i-ij^, Wi ) j ) k-,i ) j±i,k, and Wij t k-,i,j,k±i- Because the 
spins take continuous value, the summation for spin value turns into the integration 



/oo 
/ {a) da; 
-oo 



and then the normalization factor 



Qi,j,k;i+l,j,k 
+00 +00 



b 

2^ 



— OO — OO 



X exp {K [&ijk {&i s j+X,k + <J i-lJ ) k + &i,j-l,k + a iJ,k-l + + Vijk&i+ljik 

( w 2 2 



TV 



2ny b + K 



exp [if (cr iife + CTi+ij, fe) 



6 2 

x {&i+l,j+l,k + < J i+2,j,k + + CTj+lJ.fe+l + C'i+lJ.fc-l) — 2 v^y* <T *+ 1 J. fc ) 



x exp 



4 (6 + K ) 



Thus this combined term in Eqs.(f|) becomes 



&ijkWi t j^-i+l,j,k ( cr ijfc cr i+l,j,fc y °ijk&i+l,j,k) 

^ +OO +OO 

nZTl / / ^^' fc ^* + 1 'i' fc< ^' 5 "iife+ <3 "i + l,i,fc'' J ^'=+ cr i + l 

Z7T l^i j k;i+l,j,k J J 
u> wi _ OC) _ QO 

xci-jjfcexp [<j„ fc (crjj + i^ + + cr i; j_i >k + Oij^-i + ffijfi+i) + &ijk&i+i,j,k 

+°" l +l,j,fc + °"i+2,j,fc + 0"i+l,i-l,fc + Ci+lJ,fc+l + CTf+lJ.A:— ljj — 2°"ij ~~ 2 CTi+1 ' J ' 
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[K (°"j,j+l,fe + a i-l,j,k + a i,j-l,k + &i,j,k-l + a i,j,k+l + a ijk + 



2 (b + K) 

—&i+l,j+l,k - a i+2,j,k - &i+l,j-l,k - &i+l,j,k+l - &i+l,j,k-l) + b ((Tijk + CTi+l,j,k)] ■ 
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